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PREFACE 



A number of thermodynamic charts for internal combustion 
processes have been prepared in the last twenty years. None 
of those presently available, however, cover the pressure, 
temperature and fuel-air ratios of diesel cycles. 

The purpose of this project is the preparation of one 
thermodynamic chart of a typical diesel fuel-air ratio, 
covering the usual temperature and pressure ranges encountered 
in the combustion portion of the diesel cycle. The work was 
undertaken at the suggestion of the Internal Combustion 
Engine Division of the Bureau of Ships, which furnished 
the author with the necessary data on pressure and temper- 
ature range, and fuel-air ratio. 

During the period from November, 1950 through May 1951, 
the author made the analysis and computations necessary for 
the project and constructed the thermodynamic chart. Inves- 
tigation was further made toward a simplified method of 
attack for the problem of construction of the low pressure- 
low temperature region of this chart for determination of 
the thermodynamic properties of the combustion products in 
the remainder of the diesel cycle, i.e., the adiabatic 
reversible expansion, and the opening of the exhaust valve 
and release of pressure to atmospheric. 

The author gratefully acknowledges the timely assistance, 
advice and guidance of Professor Gilbert F. Kinney of the 
Department of Chemistry and Metallurgy of the U. S. Naval 
Postgraduate School throughout the entire project. Acknow- 
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ledgements are also extended to Professor Paul J. Kiefer 
of the Department of Mechanical Engineering of the U. S. 
Naval Postgraduate School for valuable suggestions and 
advice anent the project; and to Lieutenant E. W. Laessle, 
USN, for assistance in drafting the charts. 

The entire project was conducted at the United States 
Naval Postgraduate School. 
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CHAPTER I 



INTRODUCTION 

The large fraction of combustion processes of interest 
to the engineer involves burning of a hydrocarbon of fairly 
moderate volatility in air. For purposes of analysis by 
thermodynamic chart, an internal combustion engine cycle 
may be divided into (A) low temperature regions where 
chemical dissociation does not occur, or negligible in 
importance, and (B) high temperature regions, reached in 
combustion of the fuel, where chemical dissociation may 
become of appreciable magnitude. 

Charts prepared by H. C. Hottel, G. C. Williams and 
C. N. Satterfield in 1949 (2) , cover both the above regions 
for the gasoline engine cycle, and are satisfactory use 
with diesel engine cycles in region (A). They do not, 
however, cover the temperature range and pressure range of 
the diesel during combustion, i.e., in region (B). 

The thermodynamic chart presented in this project 
covers the pressure and temperature range of the combustion 
portion of the diesel cycle for one typical diesel fuel-air 
ratio. It has been constructed on the same material and 
temperature bases as charts for gasoline engines in region 
(B), and thus may be used in conjunction with charts rep- 
resenting the unburned mixture, as prepared by Hottel, 
Williams and Satterfield (2) . Values of equilibrium con- 
stants and thermodynamic property values used here are 
consistent with those used by Hottel, Williams and 
Satterfield (2) . 
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CHAPTER II 



BASIC CONSIDERATIONS AND DATA 

1. Data Furnished by the Bureau of Ships 

It was suggested by the Internal Combustion Engine 
Division of the Bureau of Ships that the thermodynamic 
chart be constructed for the following: 

Fuel-Air Ratio = 0.5 (200% theoretical air). 

Pressure Range: 800 p.s.i. to 1500 p.s.i. 

Temperature Range: 3000 ° Rankine to 1 + 500 ° Rankine. 

2. Chemical Dissociation 

At the high temperatures reached in the combustion 
process of the diesel engine there may be many different 
molecular species present in varying amounts, due to 
chemical dissociation of the working fluid. The requirement 
of chemical equilibrium necessitates consideration of the 
effect of the species on the thermodynamic properties of 
the burned mixture. There are ten molecular species (N 2 , 

O 2 , CO 2 , H 2 O, NO, OH, CO, 0, H 2 , and H) present in suf- 
ficient quantity to require consideration in this develop- 
ment. In the course of the project, it was found that the 
effect of the minor species (NO, OH, CO, 0, H 2 , H) was of 
definite importance in the temperature and pressure range 
involved. In Figure 1 is shown the effect of temperature 
on the relative magnitude of these species for one repre- 
sentative pressure. It is apparent that the magnitudes of 
these minor species become appreciable at the higher temper- 
atures, and cannot be neglected with propriety. Particularly 
is this true of the species NO and OH, which frequently are 
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neglected in determinations of burned mixture composition. 

3. Ideal Gas Law 

Computations involving composition of the burned 
mixture together with critical pressures and temperatures 
of the component species indicate that the burned mixture 
behaves as if its critical pressure and temperature are 
about 30.3 atmospheres and 320° Rankine, respectively. At 
the elevated temperatures of the chart deviations from the 
ideal gas law are insignificant. Therefore, the chart is 
based on the equation of state for ideal gases. 

4. Fundamental Data 

The fundamental data required for calculation of this 
chart are: (a) equilibrium constants as functions of temp- 

erature, and (b) various heats of reaction. The author has 
accepted for his calculations the same data on internal 
energies, enthalpies, entropies, and equilibrium constants 
as did Hottel, Williams and Satterfield (2). The method of 
determining the changes in the thermodynamic properties for 
the chemical reaction in question is given in Appendix I, ' 
page l(o ; the method of obtaining the actual equilibrium 
constants used in the calculation of composition of combustion 
products in equilibrium at the various temperatures and 
pressures if described briefly in Appendix I, page /^ . 

Composition of air: on a basis of 100 moles of dry 

air, a composition of 20.99 moles of oxygen and 79.01 moles 
of nitrogen was used. In actuality, the nitrogen above has 
minor components of argon carbon dioxide, and hydrogen in- 
cluded, and might better be termed Air-Nitrogen. The above 
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composition was used in calculating thermodynamic properties; 
hut pure nitrogen was used in equilibrium calculations, which 
made but a negligible change in the amount of NO present. 

Composition of fuel: the fuel chosen was the hydro- 

carbon (CH 0 ) , which has a hydrogen-carbon ratio considered 
adequately representative of diesel fuel. 

Other constants: 

1. Gas Constant R - 1.987 cal./°K. gram-mole. 

2. Apparent Formula Mass of Air - 28.95 (which 

agrees with the above air composition.) 

3. Atomic masses from 1936 International Tables. 

4. Conversion Factor: 1 atmosphere - 14.696 p. 3 .i. 

5. Conversion Factor: 1 cal. /gram-mole - 1.80 

BTU/lb. -mole. 

6. Conversion Factor: 1 lb. -mole burned mixture - 

0.0024135 pound of air. 
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CHAPTER III 



PROCEDURE 



1. Analysis 

If definitely fixed quantities of carbon, hydrogen, 
oxygen and nitrogen are permitted to come into equilibrium, 
two of the properties of the resultant mixture may be fixed, 
and will be sufficient to completely define the composition 
of the mixture, and hence its thermodynamic properties. 
Pressure and temperature were the properties selected in 
this work. 

It was considered thRt chemical dissociation would 
produce combustion products containing the following molec- 
ular species: N 2 , O 2 , CO 2 , H 2 O, CO, H 2 , H, OH, 0 and NO. 

The relative magnitudes of these species, at the fixed fuel- 
air ratio, depend upon the pressure and temperature. The 
method used to calculate these magnitudes at a sample 
pressure and temperature, using determined data on equilibrium 
constants, is shown in Appendix II. 

2, Material Basis 

The material basis of the chart is the quantity of 
material which contains nitrogen and oxygen equivalent to 
one pound of air. 

For complete combustion: 

(c»X + + 3.7U 0*3 A/J = A (COj 1 - H x O +■ S.C.H- A/J 

For a fuel-air ratio of 0.5, there are three moles of oxygen 
and 11.280 moles of nitrogen, or 14.280 moles of air, for 
each CH 2 unit of fuel. The mass of air is (14.280) (28.95)= 
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413.406 pounds; the mess of fuel is (l)(l2)+- (2)(l.008)= 

14.016 pounds. The mass of fuel associated with one pound 
of air is then 0.033904 pounds. The mass basis is necessary 
for consistency because the number of moles of combustion 
products varies at elevated temperatures with variation of 
any thermodynamic property. 

3. Determination of Thermodynamic Properties 

The temperature base for the chart was selected as 
100°F. The values of internal energy and entropy, obtained 
from Table 6, Appendix A, Hottel, Williams, and Satterfield, 
(2), were converted to this temperature base as described 
and illustrated in Appendix I. To determine the internal 
energy of the combustion products at each selected tempera- 
ture and pressure, the procedure is to multiply the mole 
fraction of each molecular species in the products by its 
corresponding value of internal energy, corrected to proper 
temperature base, sum up these products for all ten species, 
and convert to the proper material basis for the chart; thus 

E . = US: ■ 

Entropy involves not only an additive factor similar to that 
shown above for internal energy, but also the entropy of 
mixing; that is, the entropy of the mixture is 

= ^ Si Si Pi) 

where ^ is the mole fraction of each species, 5-. is the 
molar entropy of each species at one atmosphere pressure 
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and each temperature, corrected to the proper temperature 
base, R is the gas constant, and P; is the -partial pressure 
of the components, in atmospheres, and also equals P } 
where P is the total pressure of the mixture in atmospheres. 
The entropy is corrected to the proper material base 

for the chart. The enthalpy of the mixture is obtained from 
the thermodynamic relation 



The volume at each pressure and temperature is obtained 
from the ideal gas lav/: 



where is the total moles of the mixture at each temper- 

ature and pressure, R' is the gas constant (1545 BTU/°R 
lb-mole) and -f- is the total pressure in pounds per square 
inch. This volume is then converted to the proper material 
basis of the chart. 

4. Description of Chart. 

The thermodynamic properties represented on the chart 
are as follows: 

(a) Temperature -T, degrees Rankine. Nearly horizontal 
solid lines at 200 degrees intervals. 

(b) Pressure - -p- , pounds per square inch. Diagonal solid 
lines at 100 p.s.i. intervals. 

(c) Volume - V, cubic feet per pound of original air. 
Diagonal dashed lines with greater slope than the pressure 
lines, at 0.1 or 0.2 cubic foot intervals. 
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(d) internal Energy - E, BTU per pound of original air. 
Ordinate of chart. 

(e) Entropy - S, BTU / '°Hankine per pound of original air. 
Abscissa of chart. 

(f) Enthalpy - H, BTU per pound of original air. Nearly 
horizontal dashed lines, superimposed on the chart to elimi- 
nate the necessity of using the relation H = E + RT. 

5. Net hod of Bead n^ the Chart. 

To illustrate reading the chart let it he desired to 
determine the thermodynamic properties of the products of 
combustion of one pound of air with (CHg)^ vapor at T = 

3600° Rank Ine and p = 1200 p.s.i. The original mixture 

• mT!.- 

contains 100* excess air. Locate the intersection of the 
lines corresponding to the above tr-m'oera ture and pressure, 

i V 

and read direc + l that the volume is 1.15 cToic feet; the 

az < O: 

internal energy is 711.8 BTU; the enthalpy is 953.2 BTU; and 

14 

the entropy is 0.285 3TU/° Rankine. Interpolation between 
temperature and enthalpy lines is substantially linear. In 
actual cycle analysis employing the chart, the known ther- 
modynamic properties v; ill be more probably be pressure and 
either internal energy or enthalpy, but the procedure in 
read’nn the chart remains the same. 

6. Use of Chart. 

The chart is designed for use in conjunction with diagrams 
represen 4 at ive o" the unburned mixture existing prior to 
combustion, in analysing diesel cycles. Such diagrams have 
been prepared by TT ot 4 el, 7/ill ians , and Satterfield, (2). The 
author’s chart indicates the tw^r mo dynamic properties of the 
burned mixture as a result of the combustion process. 
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CHAPTER IV 



A DISCUSSION OF THE SYSTEM BELOW 3000° RANKINE 

1. Basic Requirements. 

To permit determination of the thermodynamic properties 
of the burned mixture as a result of isentropic expansion to 
a point where the exhaust valve opens, a low pressure - low 
temperature chart, of substantially the same entropy limits 
as the author’s chart, is required. The approximate temper- 
ature and pressure ranges for this chart are 900° Rankine 
to 2900° Rankine, and 14.7 p.s.i. to 190 p.s.i. respectively. 

2. Frozen Equilibria. 

Chemical equilibrium has been assumed to exist at all 
temperatures above 3000° Rankine, where the chemical reaction 
velocity is quite high. However, when the gas mixture is 
rapidly cooled below 3000° Rankine, the composition is found 
to correspond to that at some higher temperature, at which 
the composition '’froze 1 ’ . Considering that in diesel cycles 
an isentropic path is closest to the actual path to the 
system below 3000° Rankine, the pressure effect on composition 
was taken into account by the assumption that the composition 
of a mixture below 3000 ° Rankine was that of the mixture at 
3000° Rankine, at the same entropy. Experimental validation 
of this assumption is quoted in Chapter V, Robinson ( 5 ), 
and is consistent, in this temperature and pressure range, 
with the value of 2880° Rankine, used by Hottel, Williams, 
and Satterfield (2). 
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3. Determination of Thermodynamic Properties. 

The composition of the burned mixture at 3000° Rankine 
and at various pressures is determined as shown in Appendix 
B. The entropy of the mixture at 3000° Rankine and the above 
pressure is determined as described in paragraph 3 of Chapter 
III. As stated in paragraph 1 above, the composition of the 
burned mixture at any temperature below 3000° Rankine is that 
of the mixture at 3000° Rankine at the same entropy. Thus, 
selection of any temperature T in this region fixes two ther- 
modynamic properties - temperature and entropy - and these 
are sufficient to define the remaining thermodynamic proper- 
ties of the mixture. The internal energy and enthalpy at 
each temperature T is then obtained in the manner described 
in paragraph 3 > Chapter III, using the mixture composition at 
3000° Rankine at each entropy. To obtain pressure (p) and 
volume (V) , the equations for the isentropic process by which 
this region is entered in the cycle, from the region above 
3000° Rankine, are used. These are, at any temperature T and 
entropy S: 



where p^, and V^ are the pressure and volume at T^ = 3000° 
Rankine, and k is the mean specific heat ratio between the 
temperatures T and T^. Values of the specific heat ratio 
for the (CH 2 ) x fuel with 2005 b air, at any temperature between 
300° Rankine and 4000° Rankine, may be found in Table 8, 

Gas Tables (3)* The method of obtaining the mean specific 
heat ratio between any two temperatures, when the specific 
heat ratios at these limiting temperatures are known, is set 
forth in Kiefer, Stuart and Kinney (4). 
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The internal energy of any species at temperature T 
equals E 5 +Q V where E 5 is the heat released by cooling 
the gas at constant volume to the temperature base of the 
chart (100°F = 3H°K) , and Q v is the heat of combustion of 
the gas at constant volume at 100°F. E s is the value of E 
of the species at temperature T minus the value of E of 
the species at 100°F, each obtained from Table 6, Appendix 
A, Hottel, Williams , and Satterfield (2). at 100°F is 
the heat evolved in cooling the species, and oxygen required 
for combustion, down to 0°K, plus the heat of combustion at 
0°K (E at 0°K in Table 6 } , less the heat absorbed in heating 
the combustion products from 0°K to 3H°K (100°F). 

For the combustion of CO, the combustion equation is: 



or one mole of species CO reacts with one half mole of 
oxygen to form one mole of CO 2 as combustion products. From 
Table 6, Appendix A, Hottel, Williams, and Satterfield (2), 
the following values are obtained: 



CO + 'U o l 




For CO: E Wk = 66,769 For C0 2 : C„,.„-C 0 . K = 1737 




Then: 



E-s = = 7S435-73757 = 467 8 
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= (73757 - 66769) + + 66769 - 


1737 




r 72789 














= 4678 + 72789 


= 77467 calories/g-mole. 


The 


above is the value shown in Table 2. 




The 


value of enthalpy 


is readily obtained 


from the 


relation 


H = E + RT, Thus 


the enthalpy of CO 


at T - 3600°R 


is: 


^36>oq‘R ^-3(*00°R 


+ 1.987 * 2000 






= 7 7 H- <ol 


■+ 3974 = 81441 calories/g-mole 




TABLE ; 


III 




Molar Entropy ( at one 


atmosphere pressure 


) 


Species 


T - 3000°R 


T - 36 oo°R 


T = 4500°R 


o 

o 


20.166 


22.522 


25.791 


h 2 o 


15.910 


17.875 


20.68 


Air-N 2 


12.754 


14.148 


16.072 


02 


13.490 


14.956 


16.993 


CO 


12.924 


14.336 


16.279 


H 2 


12.061 


13.537 


15.397 


H 


8.34 


9.26 


10.36 


OH 


12.40 


13.74 


15.62 


NO 


13.21 


14.65 


16.61 


0 


8.34 


9.26 


10.36 



The units of the above are calories/g-mole ° K , or BTU/ib- 
mole °R. 
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The values shown in the preceding table (Table III) were 
obtained from those in Table 6, Appendix A, Hot tel, Williams, 
and Satterfield (2), merely by subtracting the entropy value 
at the temperature base of the chart (100°F = 311°K) from 
the entropy corresponding to the temperature in question. 

TABLE IV 

Thermodynamic Properties of Equilibrium Mixture of 
Combustion Products From Which Chart Was Constructed: 



T 


P 


E 


H 


S 


V 


°ft 


pst 


BTU 


BTU 


BTU/°a 


cu.ft 


3000 


800 


548.376 


761.731 


0.25860 


1.4399 


3000 


1200 


548.356 


761.711 


0.22988 


0.9599 


3000 


1500 


548.347 


761.701 


0.21385 


0.7680 


3600 


800 


711.951 


953.289 


0.31384 


1.7281 


3600 


1000 


711.855 


953.191 


0.29724 


1.3825 


3600 


1200 


711.777 


953.110 


0.28499 


1.1521 


3600 


H00 


711.720 


953.052 


0.27403 


0.9875 


3600 


1500 


711.687 


953.017 


0.26909 


0.9216 


4500 


800 


976.566 


1297.191 


0.39238 


2.1633 


4500 


1200 


974.417 


1294.959 


0.36334 


1.4422 


4500 


1500 


973.242 


1293.749 


0.34733 


1.1536 



The material basis consistent with the above values 
is one pound of dry air (1/28.95 pound-mole), plus the 
corresponding fuel. 
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Method Of Calculating Composition Of Combustion Products 

In Equilibrium 

The large number of molecular species present in the 
combustion equilibrium mixture makes it necessary to organize 
the method of calculation to minimize trial-and-error . 

There are four kinds of atoms present (C, H, 0, N) . 

These are distributed in ten species (COg, ^0, O 2 , , CO, 

H 2 , H, OH, NO, 0). Available in relating these species are 
the six chemical equilibrium equations shown in Table I, 



Appendix I. They are listed again here for convenient 

reference, with the number of moles of each species being 

/ 

represented by its chemical symbol in parenthesis: 


1 


K, = 1 




t v 2 


K^WPh,!* 




E v 3 


k 3 = <U,/ c = 










E v 5 


K, -• Vf £ 


(H) /( H y> *if% 




= p °" K\/ 

' • w 2 o 





In the above, P is the total pressure in atmospheres, 
is the partial pressure in atmospheres of each molecular 
species, and the symbol £ represents the sum of the moles of 
all species present. 



( 20 ) 



As shown in Chapter III ( 4 ), with a fuel-air ratio of 
0.5, or 200% air, there are 3 moles of 0 2 and 3.76(3) 

11.280 moles of N 2 present, giving 1A.280 moles of air per 
mole of fuel. If one mole of carbon and of hydrogen in the 
fuel is chosen as a basis, there are available the following 
four equations satisfying the material balance of the four 
atoms present (C, H, 0, N) . These equations are, on a mole 
basis, at equilibrium: 

EL<i. 7 (co z ) ■+ (co) -- /. ooo 
E a . 8. (n z o) + (HJ + '/z/ThWohVJ = / ooo 
£ ^ 8 . (coj 4 - (0,3 + K[(^o)+(cd) +(ou) + (no) + ro)j = 3. o oo 
10. (/O +^(HO)~ll.2QO 

The equations 1 to 10 inclusive give sufficient relations 
to solve for the composition of the ten molecular species in 
the combustion products. 

Procedure: The first step is to divide the ten species 

into major components (C0 2 , H 2 0 , 0 2 , N 2 , CO and H 2 ) , and 
minor components (H, 0, NO, and OH) and to consider for the 
first approximation of composition that only major components 
are present. Then equations 7 to 10 inclusive become 



% 7A. 


(co,) + (co) * /■ ooo 


8 A. 


(H z 0) ■+ (hJ -/.O oo 


Ecj. 8 a. 


(coO 4- (oO 4 '4 {(K6) + (co)J = 3.0 00 
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E.^. IOA (M,) = 280 and 

E^. II. £ * (co x ) + (h,o) + Co,) +(nJ + Cco) + Cm,) 



By judicious use of the above five equations together 
with equations 1 and 3, a first approximation of the number 
of moles of the major species present may be obtained. Then, 
using these results, plus equations 2, 4 , 5 ana 6, a first 
approximation of the minor species is made. The values of 
the minor species so obtained are put into the original four 
equations 7, 8, 9, 10 and these equations in turn solved in 
terms of the major species. A second approximation is now 
carried out, in precisely the same manner as the first 
approximation, to give corrected values of the major species 
based upon the presence of minor species. Again values of 
the minor species are determined with the aid of equations 
2, 4, 5 and 6. If the results for all species now satisfy 
equations 7, 8, 9, and 10 the computation is completed; if 
not, further approximations are carried through until these 
equations are perfectly satisfied. 

To illustrate the method, a sample calculation will be 

t 

performed, using representative values of temperature and 
pressure. 



Let T- 3000°R Let -73* 800 P . S ., , P= 54.4 atmospheres 
From Table I, Appendix I, at T= 3000°R: 



K, = 3.406 

= 9.165 io “ 8 

K 3 = 1.816 10" 10 



K, = 3.981 10" 10 
K, = 1.014 10-*- 

K 4 = 5.053 10- 6 
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(cc^GO 
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toOCwJ 

<w z o) 2 
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CM 



wo 



or 



* '.^7 '0 S [/^] 




Adding equations 7A, 8A and 10A: 

Ea.12. (coj+Ccdt+C^di+CnJ+M = I3.zqo <i*d 
E<*. '3 £ = I3.Z8Q + <&) 



Multiplying equation by two and subtracting equation 7A: 



t<j M. (COj + 


Ch^o) + 


iCoO = 5.000 






From equation 10A, 


w = 


II. 280 






1ST 


APPROXIMATION 










.L 


A: 








/.500£>l 


1. 5000 12 7 


(GO 


£=15.280 + (o^) 


& '3] 


/ ^.7804,1 


1+ 780012 7 




PC 

VFCOJ 




0.428 


0.4-28 




347/0^ 


HA 3 1 


0.515 to' 5 


0.515 /0~ s 




60 




0.00 000575 


0.00000575 


<K) 


(h 2 oW.ooo-GO 


El- 8 a] 


0.44994425 


0.444444 25 


6tp) 


(COj) » 5.000 - 2 (o>(H 2 0) 


Ch <+) 


0.44 8 78575* 


0.94 94804 


coo 


(co) = /.ooo - (CoO 


fil-TA) 


O.OOUIM25 


O.OOOOI96 


fco) 




B*.iJ 


• 


3.doo 
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* 



£<2.5. (h) = CM -0O Vi Vf 

~ (1014 IO' f )( 0000057S) V% ( i±J_&OOU]'\ / - 



£0.6. 



U.2 



= 0.00000 1266 



(°«)* ( ^ • ggfc-Jf 

-(5:053 • . / ^.7800)37 ^ 

(o.o000057j)*^ '• 5'^.4 ' 



~ 0.001100 



Ovo) = 00 - • 00 v * - v/X 

ch;> » p 



- (9. / 6 5 • /o' ? ) /°.<»^W2i) 300 ^ f/^ ysootzi W' 

( O . OOOOOSTS ) ^ ^ 



-- 0.0280 



E ,.4 (Ol = <K»} ■ f i 2 } • i, 

^ (MO F 

- (3 c}Q| . ,0-'°)- f '^!-L eoOU, N ) 

(0.0000 057.5) ^ / 



- O.OOOOI8 72 

Substituting these values for the minor species in 
equations 7, 8, 9, 10: 



£^7 (coj -t- (co) - I.000 

L^. 8 Ch 4 6 ^ + (H*) - / ooo - ‘4 +<oh>] 

= 0.999 99 9 37 

EL^. ^ (coj +(0^ + ^[(H i O^(co '0 = 3.000 - 4 E oh) +Ca/o 3 4 (o)J 

~ a.i8 I 
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E. % . 10 (N t ) ^ 11.280 - '/z (NO) 

= 11.266 

Now, following the same procedure as before: 
Adding Eq. 7, 8, and 10; 



UA. CCOj) t (CO) ->• (H z O) -+(H Z ) + (N 2 ) = /3.2.65'44 c |4 and 

E t /3a $ - /3.1654494 + (0,) 

Multiplying Iq. 9 by two and subtracting Eq. 7: 

EL^. /‘/A (co z ) -*■ (up) + a(oj = 70 882 



2ND APPROXIMATION 



(0.) 




/.moot 


1. 48S727 


(oJ 


£=/3.X6s<m4+<b;> 




1 4. 155HC, 


/4. 75-/ / 7 8 




Jl 

Ykw 




0.427 


0.428 




W /«.o)'«47 v’JSj 


[Ei 5] 


o.5is/o' s 


0.5 75 ■/6 s ' 




(H.) 




0 0000057.5 


O.OCOOQ57S 


(H,) 


(8,0) = 0 444 4 4 4 3 7 - (H,) 


CM] 


0.74744362 


0 99944362 


£4,0) 


(CQ>4770882- 2ftO-(H 4 0) 


El MA] 


0 7714184 


0.4447804 


(CO,) 


(co) = /, oo o -Ccoq 


fay 


0-0085816 ‘ 


0. OOOOI46 


(co) 


(C0)(H,0)/ _ 3 ^ 06 

/(coj60 - 


fay 




3 406 
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Again determining the minor species from Eqs . 5, 6, 2 and l+: 



oo. O.OOQOOIJ66 

(oh) = o oonoo 

(NO) = 0-03.8 oo 
(o) = 0.0000/872 

The values of these species have not changed from those 
obtained in the first approximation, and hence all values 
satisfy equations 7, 8, 9, and 10; if they had not, a third 
approximation would have been required, performed in the 
same manner as the above. 

Final Results: 



fcoj 


=-0.9^99804 


(W) 


- 0.999 44342 


(o z ) 


= A 48 5721 


CaO 


= J/. 266 


(CO) 


^ 0.0000 G 


0-0 


= 0.00000575 


Cm) 


~ O.OOOOOI246 


( oh ) 


- 0.001100 


(no) 


= 0.02800 


Co) 


- 0.00001872. 


A 
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